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SUMMARY 


A photovoltaic power system with a * 

A dual battery current centre! eencept < * pabm *>' ,a >na ^^- 

te either supply or acce ^ power P P " >poaed - whlch '"»>>les the battery 

conditions. A simulation of the power systeiTVTT e " vironment and load 

control, is developed and evaluated The g the ba * tel 5 r current 

of the battery control concept to switch the' bLt'” lMO f demonstrated viability 

mode and back as required by load and » ’ r ° m a clwr 8c to discharge 

system operation Is demonstrated over t hI en^T‘ C ° nditl0n8 - An “ooeptable 

system sensitivity, bandwidth, and damninv I S ° la ™nge. Additionally, 

trol are shown to be acceptable for fo^T^ racteristios <* »o tottery con- 

P we for a projected hardware implementation. 




to more usable energy forms. One mettled is the dl " 0t solar 

eolation to a do vol.age «, current by the use of ~ nversl «" »f solar In- 

then be ’'inverted" to ac power comnaHM <*v 1 CeU8, 11118 dc can 

The Terrestial Photovoltaic Project Office ™ ri ° US ® C machines loads, 
veloping a 10 kVA Inverter/Contm.^ K Pflrt of a 00 E ^ort, Is de- 

Power System (PPS) at the Lewis R es *ea U h *** ^ Sting in 8 Photovo l tai o 
combined system (PPS) is to inv l t T (LeRC) ‘ purpose of the 

power sources to supplement and/or*** 6 ^ 6 applicatI on potential of photovoltaic 
description of the energy suppllee. A 

In a true "stand ! P °“ ^ tound reference 1. 

require energy storage capability 8 °to,' ° OMinuous mode of operation, the PPS wl 

‘ode of low solar Insolation and accept exe^vaTlW 0 " 1 " T ^ Pe 
during peals sunlight hours in order to ontimi . ^ * P ° Wer fr ° m the SoIar arr * 

Power capability. Such storage capability is typicall^^ eXP * n8iVe SoIar arra y 
Simplistic batterv storage systems h«v h ly aecompU8hed with a battery 

tions by simply shunting^^ d ~ power PPS install! 

Optimization of these designs for lLe!h« V ^ ° UtpUt termlnal «- 

temperature is difficult. Also the butt ^ " S ° Ur insolatlon and operating 
ta over change and charge depletion. re^n ^s a 10 PO ’ S ‘ b “’ ““ 

Battery Current Control (BCC) eonceot aa intludes " m “re sophisticated 

native to this simplistic design. ‘ to ‘" PPS W™- “ «" alter 
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0) an wrTof^lJoX’ “) tteTo “* m>i ° r bl °° kS: 

power ey8lero controller (4) ? ** 

trol (DCC), and ( 8 , the battery current control and batterv C ° n ' 

reetly convene eolar Insolation Into dc power The .nver.^? , " TV d ‘- 

^ztzz:z:zt t mMt 

system losses. Should the oc»r !° ^ the dedicated loads and 

demand, the PSC s^L Z * 8A * to — . the 

the de bus to fr °” * he ^ 10 

exceed the demand, the PSC signals the BCC a rt * P ° W6r generfl ted b y the SA 
from the dc bus to the batterv Th . excess power is transferred 

approximate maximum power operatl^volUgTLrthTs 8 ^ 168 ““ P8 ° W “ h “ 
te in a state of full charge, a s£al is sen! m 2 DCC ^ 

dissipated through a resistance load. Functionally the BCC oTT ‘ S 
ure 2. The two dc to dc converter* oc BCC operates 88 In fig- 

from the battery. Additionally, thev act a/™” SWU ° he8 to direct the flow to and 
since at any instant of time the dc bus it P ° Wer conditi °ning or matching devices, 
tery voltage and cur rim ^ Can differ the bat- 

a m ode r o Hhe C ^ batte ry W ^d ^co nt r ert s T “ C< ~ by fo —‘*tin g 

tern mode, descrlbe’d InTe^Te IZZZZ T*** “ ^ 

PSC to handle additional control function^ mu l . , ^ made *° the 

for the analog computer simulation of the comblnedTpS T m0del “ Sed 

several experiments performed on this simulation Then * ’ reSuUs of 

to evaluate system operation for various Insolation IcTels 3t!v T* ”T ^ 
configurations. This report berlns with a ’ 1 d 1 ls ’ ^ control 

components. " ** * detWled of Urn PPS and its 


Photovoltaic Power System Description 


For the simulation the PPS was modeled as shown in fi , a „ 
bols lS given in appendix A and a summary of mortal « 3 ‘ A Ust of sym ' 

dtx B. Dynamically, the important elements are mT T e ‘ Ve " " ‘ PP ™~ 
Power inverter; (:i) the power system n f .. 1 h Ur array{ (2) the 

capacitors; (5) the battery current charge Ind d7s h' ^ ^ °'° 21 farad fUter 
battery. In this model the dUo t sen ^charge converters; and (6) the 

array open circull voltage data .node”' f “ nCt '° n ‘ 8 0bla ‘ ned rr ° m the 

variables. *»d«. model, figure a also Indicates Important system 
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3o 1 J3aZ7;:i:r;:r is r by the s ° ,ar array - *■*««•■« *«<* 

"inverts" do powerTv arvl l , , T “ d °“ rrent - ‘nverter 

constant voltage magnitude. TheTh! ° V ° ltage and current to ac power at a 
voltage and current of the d'c bus to a "do' C °" Verter " converte ” the varying do 
compatible with the battery state of charge (^C,' ** T 
meat. The discharge converter n«rf rt g d h P ° Wer transfer require- 

charge converter, only the direction ofTh * he Same function 88 the 

The power inverter and the charo- 6 ^° Wer ( current ) flow is reversed. 

rents. The input filter capacitors isolItTth™ *"**7 PUlSSting dc input cur - 
currents. Because of the storal T . ***** ^ theSe pulsatin g dc 

tor is needed at the input of the discharge ° f ^ h&tteTy ' no filter ca Paci- 

of the S A is determined by the pilot sen! ° nVerter ' ^ °P®*atiag point voltage 
able power at any ***** ** the PSC to — avail- 

detail. omponent model is now described in 


Solar Array 

current-voltage characteristic^ foHlch c^for Solar celIs - The 

and ambient temperatures arp * * for different solar insolation levels 

figure 4 . The entire ~ enta1 ^ - are presented in 

with 480 series connected ° f S ° lar celIs 

acteristic could be obtained bv a simni u fore a total solar array char- 
as shown in figure 4 The ordinate • . C S ° a ng ° f the ind ividual cell characteristic 
circuit current V,I, “ * “* la aall * d =-ort 

E SA, n- A normalization of each“by ToT T* ““ ° Pe ” C ‘ r ° U,t V0,tag *- 

current values yields a single curve represents^ l(7l[ V ° ltage “ d Short clrcu ‘ t 
tion conditions (fig. 5) 7, . f 811 tem peratures and insola- 

simulation sincl ^ * * — 


Power System Control 
U was determined (ref. i) t hm maximum 


maximum power from th ,. SA , s obta „ lud wh „„ 


b’o A 0.782:1 K 


SA, N 


* 


( 1 ) 


One control task of the PSC is to maintain the relation of (1) by controlling the 
flow of power (current) from the SA to the combined load of the battery, the dc 
load, and the inverter. However, operation of the dc load control was not 
analyzed in this study. A second control function of the PSC is the regu'ation of 
the magnitude of the inverter output voltage to 240 volts. Although the hard- 
ware to accomplish this second control function is actually part of the power in- 
verter, conceptually it will be included with the PSC in this study. The PSC is 
summarized in block diagram form in figure 6. In figure 6(a) the power inverter 
control voltage, E p> 0QNT , i s generated by a proportional-integral control that 
maintains the inverter output voltage, e A , magnitude at 240 J2 volts. The dy- 
namics of the measurement of |e A | for the feedback signal was not considered in 
this report to simplify the study. No appreciable impact on system operation by 
these dynamics is expected. In figure 6(b) the dc bus feedback voltage, E mr5! , is 
filtered by a three pole, 40 hertz, Butterworth filter to remove 120 hertz note* 
and compared to the reference voltage of (1). The error signal is filtered by a 
500 hertz, single pole, low pass filter to remove high frequency noise. The 

filtered error signal then drives a proportional-integral control to generate the 
power control voltage, E C0NT . 

Once the power control voltage has been generated, the charge and discharge 
control voltages are determined from it. The basic BCC scheme is 


■ - °' 8+0 ' 2E CONT 

AI > 0 


C, CONT So 

AI*0 

(2) 

/0.8-0.2E 

AI < 0 


D, CONT \0 

AI £ 0 

(3) 


where 


AI " l SA “ *P (4) 

Several modifications to this scheme were evaluated to determine if control im- 
provements could be made and to study the effect of various hardware limitations 

in the implemented control. These modifications will be discussed in the results 
section. 
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Power Inverter 

Trr lnverter transforms dc l«P«t power generated by the eolar array 
to single phase ac power compatible with most electrical machines Con- Y 

verston Is accomplished by approximating the desired sinusoidal Inverter outnut 
W lage by a preprogrammed seguence of do pulses. The resul.»t Ce„ceT 
Altered to remove high order harmonics. Additionally, the InverterTts L 1 
impedance matching device. The variable voltage and current oftte Z CiT 
matched to the ac load at a constant output voltage. For the m ,l,usZ m l 

witt d T"’ ‘ h ! ‘" Verter ‘ S consldered ideal power conversion device 

with dc power in eou&l to sc Dower ntif* r PH^ 4 n t . 

+ . . power out. The inverter losses and the outnut fii 

r„T er ‘r° ad 

variahi ? f , y th m ° de1, Therefore only the magnitude and phase of the 

: terror; are of it; r * toi - 

in block form in figure 3 * *"* ‘°“ i imp6dance is sh »»" 


Power Converters 

A power converter couples a dc power source to a dc power load where both 
source and load operate at differing voltage and current levels 

he charge and discharge converters are identical in operation They both 
couple the battery and the dc bus. Only the direction of current How Is differ! 
^erefore, ,h. mode, used for each device will be the same 1^“ 

tz: rssrr; c zzr ,y - r ° f th ° 

e“eTClt r P ° Wer orT * * “ ‘ a 

voltage is directly moduli^ 

e conver er are modeled by a diode and two resistors. The charge and discharge 
converters with the loss elements are shown in figure 3. ^ 
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Battery 

be th,t ihe state ° f *• - 

the current through the battery, and the bltte^^ate moTl -8 ’ 

errria: r rrr t? in ^ 7 

, yptcai lead antimony /lead calcium storage batterv f n ,wi* * 

analog computer simulation, the curves of fim lr « 7 88 y* To facihtatt 

line segments ~ t , figure 7 were approximated by straight 

ane segments. The reciprocal of the slopes and "x-axis" in a/** S 

P . . , , 6 u v against HOC. This enables batterv voltaic 

B’ to be calculated as a function of battery current I and SOC Ji ^ 
tions (Bl) to fBim Th.,o + * ll ’ A B’ 811(3 soc as in equa- 

computer to f nd E Z Z tlm0tlons are P^-mmed Into the analog 

“te ° Pr0gramm ‘ ng “ 3ing ‘ e 

rr: a ~ c h ne r :z 




b{soc} = k b i b 


ir ^ f f ei j Cen i t ° f fUH charge * r B is battery current, and K n 
t. Thus, If I B ls constant Md tte bat(ery is lnltuUy uncharg D_ 


is rating con- 
then 


SOC - Kglgt 


where t is in seconds and full charge is 100 A’ 


hr or .1.6X10° A-sec. Letting 
igt - .}, 6x10^ A-sec 


then 


% 
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3. 6X10 5 A-sec 


= 2.778X10 


r*-JL_ 

A-sec 


Simulation Results 

To evaluate the battery storage concept and to determine the best charge- 
discharge control strategy, several simulation experiments were made. These 
experiments studied PPS dynamic operation for changes in solar insolation, 

"runs" StatC ° f < * nd SySt6m 1<>ad ’ A b8SiC ® et 0f six ex P«riments called 

in til Were used 10 Perform the evaluations. The six basic runs are described 
igure 10. Two types of transients are used, the cloud cover/clear transient 

l d *V° ad *hls report . "cloud cover" translenUs sl- 

SoUr isolation from 100 mW/cm 2 (cloudless) to 6.5 mW/ 

I " 5 ( Seconds - TWs Is accomplished In the simulation by decreasing E„. 
lmearly from its 100 mW/om value to Its 6.5 mW/om 2 value. The normaUzi'ng 
or short circuit current Is then varied according to the equation 

l SA,N = 1 - 034 e sa,N ” 205.06 

A "cloud Clear" transient Is determined In exactly the same way except that Inso- 
lation is increased from 6. 5 to 100 mW/cm 2 . 

The second transient type is the load variation transient. This transient simu- 
lates an instantaneous step char e In system load impedance. Expertaento fortth 
transient types are performed at a 30 and 95% SOC to simulate a discharged and 
y charged condition. Also, the system Is assumed to operate at a constant 55° C 

tons for the's r r re! ‘ T* ^ W6re to ‘"vestlgate worst case condi- 

power "cm SA IT wT TT!®’ Cm 6 00,Klm °ns Include minimum available 
power om SA (6.5 mW/cm ) and the battery (30% SOC) and an lostaneous change 

minimum to maximum power demand (l.o to 10.0 kW load). 

There are four different groups of experiments. Each group Included ail six 
run c.ndu 1 Since ^re „ . Urge number o, experiment, the stouUtto ro- 
be wUl ngUre " aOCOrd ‘" S 10 *"<• flsnre num- 

ber This will facilitate comparison between similar transients. Additional com- 
ments are included to Indicate modifications made to the BCC. Note that not 
every run in every group la presented. For example, results corresponding to run 3 
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;;:;°r a T d for “ y group - Those rosuits ' ™> other om i, ted 

tests, although completed are not included here for two reasons Fir«t tKi 

findings. CaSe> Wel ' behaved Md do d “-"Onstrate any Important or unuanal 

■mi .f * luT ^ ° f ^ W “ con,plo “ d fo r «>o BCC of equations (2) to 141 
This Is the basic charge control system and this set of runs Is called group , 

periments "nTs ^ reCOrd<id ^ “ Per “" <i "‘ SyS,em 'or «U ex- 

TheSe Slx rjns s “^od as a baseline for the remaining data All 
recorded traces indicated that the total system is weli behaved I., 

ionts are seen in the control signals and in run E T. I 

limits when tlie steps in system load . *’ Ep - °° N T’ exceeds design 

14 summarize how well ZT ““countered (fig. 12) . Figures 13 and 

s : scr sr 

BCC changes the battarv f °" S from the operatln S line occur when the 

Much of this digression is directly related tath'T’t '“"T ^ dl!, ° hargG - 

PSC Fi Jiiro i •} y e atet * *** integral control action of the 

• ifc>urc 13 corresponds to runs l anH 9 u4ifi r* 

4 and 5. 2 “ hlle flgure 14 corresponds to runs 

The effect of the two low pass filters that follow the BCC s. m 
gated. Nominally, as in the aroun i a BnB - . , , ' 6 > wus invest!- 

pole with a 1.0 hertz ' ‘° W P “ SS fUters werG sl ” gl “ 

to 0.5 hertz filters “eb“ u a et o ^ P “ “*•" Gba “ g “' 
line characteristics are given in figure 15 f rU ” S simulated - The operatin 
end 5. In all cases the group i clm/! , ^ 1 “" d 2 “" d «*“«> 16 for runs 

control action than group 2 The a T f 0 " m ° r ° Cl0Sely “PProximates the idea 

lents than that of“Z! Thu is .t” 1 *"? * ^ ‘°" ger bUt 

wioth and additional dlping ™ s ^^,1 ‘ CW ° r COnlr °' 8y8tCm ba “ d - 

bility and the resultant contL ‘ S bG " etl ° l81 40 8y8 ‘ em Sla ' 

of the 0. 5 hertz characteristic is includ ce P able - Thus the implementatio 

r .. included in the studies of trrouo ^ rpp i , 

oi equations <2, („ ( 4) requirea measurement of I and I « .! ° n r °‘ 

can be large dc currents mono. a , lsA a ° d J P‘ Slnce these both 

ft currents, measurement limitations in j 

peeled. The transients of groun -I s„,d„ ,u hardware arc to be ex- 

and a tf,% measurement error tvni »l * St ’ Parate effects of an active hysteresi 
Of M of equation (4). The hvateT i , “ 0MUrernent cffect «. on the measurement 
.neasurenen: error modification to cqumUn^tu 81 '" ShOW " ‘ n th<- 
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AI- (1.0t0.05)(Ig A - I p ) 


(5) 


c 

*• 


(or ITT 18 “ d 19 " h0W U " ! Uf(ecte 0f 0,0 ^“teresls O" control performance 

ZI Z ToTJT ° ,ear ' r ° 8 ™>- E “ h "■“» compares » hys- 
eresis (w - 0.0 A) to hysteresis widths w -- 2.G A and w - 2.5 A Note th 

for a cloud clear transient with w - o 5 A th 0 nrr u 

ivpr DD n , 0 Af the Bcc voltage saturated and control 

.... ! WaS Ch In8tabilit y could result In an unacceptable SA 

condition and a system shutdown and must be avoided 

on J‘pp 7 Lf I T W the effettS ° f the -"“““-meet error of equation (5, 
the PPS Again cloud cover and cloud clear transients are shown that compare 

co trol performance for no error and ±5% measurement error operation Ftg- 

ITZ “ ,OS3 t ° f r Stem °° ntrol f0F « would he unacceptable 

the Tl f , problems highlighted In the group 3 transients Indicated ' 

the need for an aUernate BCC control strategy. An alternative, or modlfiedtn- 

onheMcT ?“J", Sln “"‘ r 10 the bM ‘° 00C ‘ roi - the m “*ncd version 

now E TZ t T"™ *"* Ch “ 808 the bMter ^ of operation ,s 

CONT r 811 ^1. The modified BCC In equation form is 


'C, CONT 


"D, CONT 


^• 8 + 0 - 2E CONT 
1 0 


E, 


0.8 - 0.2 E 


CONT 


'CONT > 0 
E CONT x 0 

E CONT < 0 
E CONT * 0 


( 6 ) 


C7) 


The bas c se of runs was completed for the PPS with the modified BCC. This set 

forT, Br ° UP 4 - ' n,e BCC “ ,ter b «*point waoforlsttos 

for this group are again 1.0 hertz. Overall system control war .oceptable for all 

" Pt "“" 3 1 and 2 - From **>« simulated responses. It Is apparent that more 
damping would be required In the actual system. The runs 1 and 2 operating ,T« 

:rr rl r uo : ot group 4 <ng - 22> sh ° w - ioss * «•*•* »<> *^u c*i lnt » lo ; 

use ation levels. Runs 4 and 5 (flg. 23, show some harsh transients" butoveraU 

2 how 0per * ° n 8 acceptable ‘ 11,0 control performance Indicated In runs 1 and 
2 however, ts not acceptable even though the ,„.s of control occurs at a low ,n 

Z ftKUre 24 (B,t ’ hlat ° r V of run 1 , the sustained oscilla- 

tion of several system voltes and currents. It was discovered that the r.ngl of 
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E D, OONT was 1116 Smiting factor in the above control problem 
tlon of the BCC to 


A slight modi nea- 


t’ 


10.8 + 0.2 Econt 
C, OONT _ \ 0 


E CONT > 0 
E CONT i0 


E D, OONT s: 


fo.75 - 0.2 

1 ° 


e oont < 0 
e oont * 0 


( 8 ) 


( 9 ) 


which lowers the range of the discharge converter voltage eliminated the loss of 
control at low insolation levels as is demonstrated in figure 25 which is also in- 
cluded in group 4. Some cycling remains, but it damps out in less than 2 seconds. 

However, overall system control is maintained even under a low insolation eondi- 
tion. 

With the updated BCC no current measurements are required. Thus, the hys- 
teresis and measurement error problems are avoided. Additionally, since 
is generated internally to the PSC, a good signal will ^ available for control . °? n NT 

spite of this, a cloud cover transient was run with a 20% deadband on E^™, 

(ilO c centered on E^^ = 0.0). The deadband, which is much larger than would 

be anticipated in practice, has a small effect on system performance as is seen in 
figure 26. 


CONCLUSIONS 

A simulation of a PPS with battery storage was developed. An analysis using 
this simulation demonstrated the viability of a BCC concept that effectively switches 
e attery from a charge to a discharge mode as required by system load and en- 
vironmental conditions. A modification of the BCC concept to incorporate voltage 
rather than current feedback was shown to acceptably control the PPS at all but 
very low insolation levels. A second modification to the BCC which lowered the 
voltage range of the discharge converter enabled the PPS to operate at even these 
low insolation levels. Additionally, it was demonstrated that by reducing control 
system bandwidth via the BCC filters, that more overall system damping could be 
obtained. Hnally, it was shown that the control implementation was reasonablv in- 
sensitive to limitations in the expected hardware control implementation. 
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appendix a 

SYMBOLS 


D 

derivative operator, sec”* 

Df } 

derivative operator, sec”* 

E 

dc voltage, V 

e 

ac voltage, V (60 Hz) 

f() 

system function 

I 

dc current, A 

i 

ac current, A (60 Hz) 

K 

battery rating constat, %/A 

M 

reciprocal slope 

SOC 

state of charge, % 

V 

resistor voltage, V 

w 

hysteresis width, A 

Z 

ac load impedance, ohms 

e 

load impedance phase angle, ] 

1 1 

magnitude 

Subscripts: 

A 

inverter output 

B 

battery 

BUS 

dc bus 

C 

charge 

CONT 

control 

D 

discharge 

ERR 

error 

F 

filtered 

IN 

hysteresis input 


1ST 

intercept 

L 

load 

N 

normalizing 

OUT 

hysteresis output 

I» 

power inverter 

UK K 

reference 

SA 

solar array 

SKI 

set point 

1 

converter output 
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appendix b 

MODEL EQUATIONS 


*B " *B, C " r B, D 


E B = ^B^B’ *> c > 
D{ SOC } = Kglg 


M C = f M, C* 900 ) 

(B4) 

M D = f M, D^° C ) 

(B5) 

E B, INT. C 15 f INT, C<?° C ) 

(B6) 

E B, INT, D = f INT, D* 900 ) 

(B7) 

E B, C = E B, INT, C + M C*B 

(B8) 

E B, D = E B, INT, D + M D r B 

(B9) 

e b =/ E b.c* I B><> 
1 E B, D* r B < 0 

(BIO) 



Battery Discharge Converter: 
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E «£F 0 0125 L SA,N 


d < e p. cont. ,» ■ ■ Ep ' ooNT ;,!r ,coNT,F 

0 a l0«7 


E P, CONT = E ERR ' v/ Q "ERR 1 


+£' e eri 


dt 


ERR “ | C A|SET " e A 


Battery Current Control: 


D! E c CONT I - f&PONT I E C, PONT, F 
C.CONT.F «. 169 


E, 


D „1 = D, PONT a D. CONT. F 

1 D, CONT, F' ”~7 *~ 

0.159 


C, CONT 


{ 0.8 + 0 . 

. 0 


2 Emxrrr X > 0 


J CONT 


XiO 


E 


D, CONT 


I Y - 0.2 E^ONT X<0 
0 XiO 


Basic BCC Version: 


X - AI - Ig^ - Ip 
Y = 0.8 


Updated BCC Version; 


X E CONT 
Y -0.8 


(Mill) 

(B32) 

(B33) 

(B34) 

(B35) 

(B36) 

(B37) 

(B38) 

(B39) 

(B40) 
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Figure 10. - Basic experiment description. 


Croup 

no. 

Run 

no. 

Figuri 

no. 

Figure 

tyw' 

Modifications 

1 

6 

1*2 

4*5 

12 

13 

14 

T. T. 
O.L 
Q,< 

None (task BCCi 
None (task ICCi 
None .task ICCi 

2 

1*2 

4*5 

15 

16 

O,* 

0. L 

0.5 He BCC litter cherpcterlstics 
0. 5 Hz BCC filter f hir&riteiiiie« ! 

3 

4 

5 

4 

5 

18 

19 i 

20 
21 

0. L 
O.L 1 
1 O.L 

■ — — — Mncr i.rur#CTf r I JnCS 

O.SHe BCC; «*0. ZO, 2.5A J 

0 5 Hz BCC; u • 0 2 0 2 5 A 
1 0,5 m BCC : (1 4 0.051 Al 
0.5 m BCC; (1 * 0.051 AI 

4 

1*2 

4*5 

1 

1 

22 

23 

24 

25 

26 

O.L 
O.L 
T.T. 
T. T. 
T.T. 

Uptated BCC 

Updated bcc 

Updoled BCC 
Modified BCC 
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Figure 11. - Simuletion results summery 
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